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ABSTRACT 

Two significant progresses have been made in the past years on our understanding of hot 
accretion flows. One is that only a small fraction of accretion flow available at the outer 
boundary can finally falls onto the black hole while most of them is lost in outflow. Another 
one is that electrons may directly receive a large fraction of the viscously dissipated energy in 
the accretion flow, i.e, (5 ~ 0.1 - 0.5. The radiative efficiency of hot accretion flow when these 
two progresses are taken into account has not been systematically studied and is the subject 
of the present paper We consider two regimes of hot accretion model. One is the advection 
dominated accretion flows (ADAFs) which lie on low accretion rate regime, < lOff^LEdd/c^; 
another being the luminous hot accretion flows (LHAFs) which lie above this accretion rate. 
For the latter, we assume that the accretion flow will has a two-phase structure above a certain 
accretion rate, and a simplification is adopted in our calculation of the dynamics. Our results 
indicate that the radiative efficiency of hot accretion flow increases with the accretion rate 
and is highly enhanced by the direct viscous heating to electrons compared to the previous 
case of i5 <K 1 . When the accretion rate is high, the radiative efficiency of hot accretion flow 
is comparable to that of the standard thin disk. Fitting formulae of radiative efficiency as a 
function of accretion rate for various 6 values are presented. 
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1 INTRODUCTION 

One of the most important parameters in accretion theory is the ra- 
diative efficiency. This parameter describes the significance of con- 
verting rest-mass energy into radiative energy. 



L 



(1) 



where L is the total luminosity emitted from the accretion flow and 
M is the corresponding mass accretion rate of the s ystem. Take the 
standard thin disc model jShakura & SunvaevI 19731 . hereafter SSD) 
as an example, its radiative efficiency lie s in the range 0.057 - 0.43 , 
depending on the spin of the black hole jNovikov & Thorn3ll973h . 

According to the differences of temperature and mass accre- 
tion rate of accretion flows, we now have four accretion mod- 
els which belong to two series, namely cold and hot ones. In the 
cold series, when the accretion rate is lower than the Edding- 
ton rate, M'^MEdi(= lOiedd/c^). we have the standard thin disc 
model jShakura & Sunvaevll 19731). W hen M > Mndd, the model is 
the slim disk iAbramowicz et al J 19881) . The temperature of the gas 
in these two models is roughly within the range of 10^ - lO'K. 
The accretion flows are optically thick, emitting multi-temperature 
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blackbody spectrum l lFrank. King & Rainel2002h . The radiative ef- 
ficie ncy of the former is high, and is independent of the accretion 
rate l lNovikov & ThorndfTs^ . In a slim disk, the optical depth 
is so large that photons are trapped in the accretion flow and ad- 
vected into the black hole; therefore the radiative efficiency is low er 
jAbramowicz et al.ll 19881 ; iMineshige et alj2000l ; ISadowskill200^ . 

In the hot series of model, the temperature of the accretion 
flow is almost virial. When the mass accretion rate is below a crit- 
ical value, Mc, ADAF ~ 5dl^~a-MEid with 9^ = kT^ /irigC-, we have 
the advection-dominated accretion flows (ADAFs; Ichimarulll977l ; 



Rees et al.lll982llNaravan & Yill994l;lAbramowicz et aljl995l; see 
Narayan. Mahadevan & OuataertI 1 1 9981 ; iNaravan & McClintockl 
200^ for reviewsfl In an ADAF the gas is tenuous. The Coulomb 
coupling between electrons and ions is not strong enough thus the 
flow is t wo-temperature with the ions being much hotter than the 
electrons l lNaravan& Yill995l) . The critical accretion rate Mct.adaf 
is determined by the balance between the Coulomb collision and 
viscous heating in the ions energy equation. When M <c M^r ADAF. 



' As we will state below, the mass accretion rate of hot accretion flows 
is a function of radius. In addition, the value of Mc adaf is a function of 
parameter 6. The value of Mct.adaf cited here was obtained when M(R) = 



constant and S 
present paper. 



10 . The result of other cases will be presented in the 
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most of the viscously liberated energy is stored as the gas internal 
energy and advected into the black hole rather than being trans- 
ferred from ions to electrons and radiated away; therefore, the ra- 
diative efficiency of an ADAF is very low. With the increase of 
M, more and more viscously dissipated energy will be transferred 
into electrons and radiated away until M;.r,ADAF is reached at which 
advection is no longer dominated. ADAFs have been widely ap- 
plied to the low-luminosity black hole sources including the super- 
massive black hole in our Galactic center, Sgr A*, low-luminosity 
AG Ns, and the quiescent and hard sta t es of black hole X-ray bina- 
ries (lYuan. Ouataert & Narayanll2003l : iNaravanlllOO^ : lYuanlllOOTl : 



(^Stone. Pringle & Be gelmaril 19991; llgumenshchev & Abramowicj 
Fl999i, ,200Q : ,Stone & Pringldl200lh or strong jYuan & BuI|201o[k 



iNaravan & McCUntockll2008l) . When M > M„,adaf, Coulom b col- 
lision cooling becomes stronger than the viscous heating. lYuanl 
( l200lh found that in this case, up to another critical accretion rate, 
there exist another hot accretion solution in which the sum of the 
compression work (PdV work) and viscous heating balances the 
cooling. Compared to ADAFs, this model corresponds to higher 
accretion rates and radiati ve efficiency ; thus it is called luminous 
hot accretion flow (LH AF; lYuanll200 H) . LHAFs have been invoked 
to explain the origin of hard X-ray emissions detected in luminous 
X-ray sources such as Seyfe rt galaxies and luminou s hard state 
of bla ck hole X-ray binaries jYuan & Zdziarskill2004 1 Yuan et al.l 
l2007h . 

In the present paper we focus on the radiative efficienc y 
of hot accretion flows. Using t h e dat a from lEsin et aL I l ll997h . 
iNaravan. Mahadevan & OuataerJ ( Il998l) presented in their Fig. 7 
the relationship between bolometric luminos it y and the ac cretion 
rate of ADAFs (see also lNaravan & Yilll995h . lYuanl ( 1200 ih inves- 
tigated the radiative efficiency of LHAFs and found that fLHAF is 
higher than the typical ADAF value. Both calculations are, how- 
ever, based on the "old" version of hot accretion flow models in the 
sense that it is assumed that the mass accretion rate is a constant of 
radius and the value of parameter S, which describes the fraction of 
turbulent dissipation that heat the electrons directly, is very small, 
S ~ 10"'. Both assumptions are now known to be no longer cor- 
rect after the development of the accretion flow theory in the recent 
years, as we will illustrate in detail in §2. 

In this paper we systematically revisit the efficiency of ADAFs 
and LHAFs after taking into account the new progresses of hot ac- 
cretion flow theory. This paper is organized as follows. We first give 
a brief introduction to the recent progresses on hot accretion flows 
in Section|2] We then describe our model in Section[3] Our calcu- 
lation results are presented in Section|4] The last section is devoted 
to a discussion. 



2 OUTFLOW AND VISCOUS HEATING TO ELECTRONS 



Since the original work on hot accretion flows bv INaravan & Yil 
( Il994t) , two major progresses have bee made, namely the existence 
of outflow and the importance of direct electron heating by turbu- 
lent dissipation. 

2.1 Outflow in hot accretion flows 

The evidences for the existence of outflow in hot accretion 
flows comes from both theoretical studies and observations. 
In the theoretical side, both hydrodynamical and magnetohy- 
drodynamical simulations of hot accretion flows have found 
that the mass accretion rate (or more prec isely the inflow 
rate; refer to IStone. Pringle & Begelmaiil 1 19991 for details) de- 
creases with decreasing radius, no matter the radiation is weak 



following a pow e r-law scaling of radius, M{R) oc R\ Most recently 
lYuan. Wu & Bui ( |20I2|) combined all relevant numerical simula- 
tions and found that the slopes of the radial profiles of mass ac- 
cretion rate obtained in various simulation works are very simi- 
lar, no matter the simulations include magnetic field or not, and 
what kind of magnetic field configuration and initial conditions are 
adopted. So throughou t this paper we adopt a single value for i-. 
I Yuan. Bu &Wul l l2012h investigated the origin of such a profile and 
concluded that they are because of mass loss in outflow (but see 
INaravan etai]|2012l for a different view). 

What is more exciting is that the above theoretical results have 
been confirmed by observations. One example is Sgr A*. Chandra 
observations combined with Bondi accretion theory gave the accre- 
tion rate at Bondi radius. On the other hand, radio observations put 
strong constrain on the accretion rate at the innermost region of the 
accretion flow. This rate is ~ 1 % of the accretion rate at the Bond i 
radius. Detailed modeling by lYuan. Ouataert & Naravad l l2003h 
show that M(R) oc i?" ' which is close to the ab ove-mentioned the - 
oretical results. Another example is NGC 3115 jWong etalj201ll) . 
For this source, Chandra can directly determine the density profile 
of the ADAF within the Bondi radius. Again it was found the result 
is fully consistent with the theoretical prediction. 



2.2 Direct viscous heating to electrons 

An important parameter in hot accretion flow theory is S. Early 
works on ADAFs have assumed that the viscous heating prefer- 
ably acts on the ions, i.e., 5 <c 1. However, detailed analysis 
of the microphysics in accretion flow has indicated that the elec- 
trons can re ceive a comparable fraction of viscous heating to that 
of the ions (Bisnovatvi-Kogan & Lovelace 1997; BlackmarJ I 19991 ; 
l0uataert"& Gruzinov. 1999; Sharma et al. 2007). For example, an- 
al ytical investigation of particle heating due to Alfvenic turbulence 
bv lOuataert & Gnizinovl ( iT999l) indicates that electrons can receive 
a significant fraction of the turbulent heating, provided the mag- 
netic fields are not too weak (/3 < 5 - 100, see below for the 
definit ion of /?). Similar conclusion was drawn in ISharma et al.l 
( |2007|) . who investigated the particle heating by the naturally- 
generated pressure anisotropy in collisionless plasma, the case of 
typical ADAF. Additionally, magnetic reconnection is also an im- 

S iortant source of electron heating in the turbulent accretion flow s 
Bisnovatvi-Kogan & Lovelace|[l997l ; lOuataert & Gruzinov|[l999l) . 
which can also heat the electrons significantly. Although a consen- 
sus on the value of S has not been reached, it has been generally 
accepted that i5 ~ 0.1 - 0.5. 

The above theoretical results have obtained observational sup- 
port i n the detailed modeling of Sgr A* jYuan. Ouataert & NaravanI 
l2003h . In that work it was found that to explain the observations 
6 ~ 0.5 is required. Therefore, we will focus on this value of S 
in the present work. Given the theoretical uncertainties, however, a 
large range of the value of 6 will also be considered, from 5 = 10"' 
to 0.5. 



These two progresses have major impacts on hot accretion mod- 
els. For example, with the existence of outflow, the compres- 
sion work for ions and further the efficiency will be suppressed 
dOuataert & NaravanI [199^ ; IXie & Yuarj|2008l) . Also, because of 
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the suppression of compression work, which is a crucial heating 
term in LHAF equations, the M range of Type I LHAF (see be- 
lo w for definit ion) will be smaller compared to the original result 
in lYuanI ( 1200 ih . A larger 5 obviously will increase the energy input 
to the electrons thus an increase of the radiative efficiency. 



3 HOT ACCRETION FLOWS: MODEL DESCRIPTION 
OF ADAF AND LHAF 

ADAFs exist only below Mcr,ADAF- Above Mc,-,adaf> the flo w will 
enter into the LHAF regime. LHAF is thermally unstable jYuanl 
but for accretion rate lower than another critical value, 
A^cT.LHAF, the growth timescale of the instability is longer than the 
accre tion timesca le, thus the flow can remain hot throughout the ra- 
dius ( lYuaj2003l) . Above Mcr. LHAp! the one-dimensional steady cal- 
culations show that the radiative cooling is so strong at small radii 
that the flow collapses and forms a thin disk ( Yuan 2001). Instead 
of global collapse, another possibility is that as a result of thermal 
instability, so me cold den se clumps will be formed, embedding in 
the hot phase ( lYuanl2003h . This scenario is perhaps more likely and 
will be adopted in the present paper. The LHAF solution belo w and 
above Mc lhaf are called Type I and II LHAFs, respectively dYuanl 
l200lh . We note that the id ea of two-phase accretion flow was also 
explored by other authors i Guilbert & Rees|[l988l : iFerland & ReesI 
1 19881 : iKrolikll 19981 : IWang et al.ll2012h . although none of them ever 
calculated the emitted spectrum. 



3.1 ADAF and Type I LHAF models 



The dynamical equations describing a two-temperature ADAF 
and Type I LHAF are exac tl y same, which are s umma- 
rized below /Nakamura et alj 1 19971 : iManmoto et all 1 19971 : 

'Oi, 



lYuan. Ouataert & Naiavan 2003 



dv , 

V Q^R 

dR 

- j) 



9adv 



pv 




Pi 


dp 






dR 


pv 


'ds. 


Pc 


dp 


y'dR 




dR 



2 1 dp 

-nlR -, 

pdR 

-aRP, 



(1 - l5)^vis - ^ic, 



(2) 
(3) 

(4) 

(5) 



where j is the eigenvalue: a is the viscous parameter: s presents 
the specific internal energy; ^vis(= ~'^P^^) the total viscous 
heating rate, q^^y, q\c and q^.^^ are the energy advection rate, the en- 
ergy transfer rate by Coulomb collision between ions and electrons, 
and the radiative cooling rate, respectively. The rad iation includes 
synchrotron, bremsstrahlung , and inverse Compton dNaravan & Yil 
1 19951 : 1 Marmioto et alJ[T997h . Subscripts "i" and "e" denote quan- 
tities of ions and electrons, respectively. All other non-specified 
quantities are of their usual meanings. 
The mass continuity equation is. 



M{R) = Mo 



R 

Roui 



(6) 



where M{R) = -AnR Hpv (H = c,/flK is the s cale height at /?) is the 
inflow accretion rate. lYuan. Wu & Bui ( I2OI2I) found that the profile 
of M{R) has two parts. When < 10^?^, ^ ~ 0: but outside WR„ 
s ~ 0.5. In the present paper, we set s = 0.4 throughout the flow for 
simplicity. 



In order to calculate the synchrotron radiation, the magnetic 
field strength remains to be specified. This is done through a pa- 
rameter yS = p/Pmagi where magnetic pressure p^^^ = B^/Sn and p 
is the total (gas-l-magnetic) pressure. Our calculation of syn chrotron 
and bremsstrahlung emission foUows lNaravan & Yil ( [1993) . For the 
Compton scattering part, we only adopt local Compton scattering 
processes (see Sec. [5] for discussions on global Compton scatter- 
ing). 

Once the dynamical structure (density, temperature, etc) of 
hot accretion flow is determined , we can calculate the spectrum 
dYuan. Ouataert & Naravanll2003l) . The bolometric luminosity can 
then be obtained by integrating over frequency. When outflow is 
included, the definition of radiative efficiency is a bit subtle since 
M now is a function of R. Now we define the radiative efficiency 
as. 



6 = L/M^^,c' 



(7) 



where M„et is the accretion rate at the Schwartzschild radius R^ = 
2GMbh/c^. This choice adopts the lowest accretion rate, thus gives 
the highest radiative efficiency. If we define the efficiency using the 
accretion rate at Rout, obviously it will be. 



M(/?„„,)c2 



(8) 



3.2 Two-phase accretion model (Type II LHAF) 

As explained above, the accretion flow is likely to have a two-phase 
structure when M > Mc^.lhaf (but lower than another critical accre- 
tion rate, ~ IOMct.lhaf, above which hot solutions cease to exist). 
In this type of accretion flow, cold clumps will emit optical/UV ra- 
diation, and these photons will also serve as seed photons for the 
Comptonization process, which is responsible for X-ray emission. 
The whole process is obviously very complicated, depending on 
the detailed dynamics of the two-phases accretion flow such as the 
filling factor and temperature of clumps and so on. This is beyond 
the sc ope of the present paper. Here following lYuan & Zdziarskil 
( I2OO4I) , we adopt a simplified approach, which is based on the 
Compton y-parameter. 



: 46,(1 + 4e,){T + t'I 



(9) 



where t = crTiigH is the vertical optical depth of the accretion flow. 
Although y-parameter is in general expected to be a function of 
radius, we now focus on the inner region, R < lOR^, where most 
of the radiation comes from. We use y = constant to replace the 
electron energy equation Eq. (|5j. This implies that we assume the 
gradient of y does not affect our result significantly in this region. 
Combined with other dynamical equations, we can calculate the 
dynamics of the two-phase accretion model for a given y-parameter. 

The total luminosity is calculated by a different approach. 
Hold in mind that the radiative cooling rate per unit volume q^^^ = 
F/H, where F is the radiative flux at the surface of the accretion 
flow. We determine the bolometric luminosity as. 



L = 2 J InRFdR = 4;r RH%.^idR. 



(10) 



To determine the quantity ^7,-,^, we note that at high accretion rate, 
the advection factor of electrons is nearly zero (refer to the dot- 
dashed and long-dashed curves in Fig. [3]l. Assuming that this is 
also true for a two-phase accretion flow, we set ^adv.c = in Eq. lO 
and calculate the radiative cooling rate as q^^d = Sq^is + lie ■ 
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Figure 1. The radiative efficiency of hot accretion flows defined in Eq. as a function of the net accretion rate. The open symbols represent results of ADAF 
or Type I LHAF, while the filled symbols are for two-phase accretion model. The solid (with data shown as diamonds), dot-dashed (with triangles), dashed 
(with squares), and long-dashed lines (with circles) represent 6 = 10~', 6 = W^^, 6 = OA and i5 = 0.5, respectively. The two-phase model for 6 = 10"' case is 
similar to 5 = 10"^ one, thus is not shown. The dotted curve is the radiative efficiency of standard thin disk model (6ssd = 0.1). The filled hexagons mark the 
value of Mcr.ADAF for each choice of 6. Power-law fitting results are presented in Eq. 41H and Table[T] 



4 NUMERICAL RESULTS 
4.1 Radiative efflciency 

We set the black hole mass Mbh to IOMq. For supermas- 
sive black holes, we find that the results are similar. The outer 
boundary is fixed to be Rau[ = lO^R^. So we have M„ei = 
{RiJRoMYMo = 0.16Mo. Here Mo = M(R„ui)- We adopt vari- 
ous S: S = 10"', 10"^, 0.1, and 0.5. Throughout this paper we set 
a = 0.1. Numerical simulations show that if the a viscosity is 
intrinsically the magnetic stress associated with the MHD turbu- 
lence driven by magnetorotational instability, as widely accepted, 
we us ually have a(3 = cons tant, with the constant being of order 
unity telackman et alfcOOSl) . We therefore set /? = 10. 

The results of efficiency are shown in Fig. [T] For given out- 
flow strength (s = 0.4), the critical net accretion rates (M(., adaf> 
Mcr,LHAF) are similarfor various ^: (6.3,7. l)x 10"' MEdd {S = 10"'), 
(6.2,7.1) X 10"'MEdd (6 = 10"2), (5.9,6.6) x IQ-'Mndd (S = 0.1) 



and (4.4,5.3) x 10 'Medd (S = 0.5), respectively. Several results 
can be seen from Fig.[T] 

• The radiative efficiency for ADAF is positively correlated with 
the mass accretion rate, as expected. When M„et < 2x 10"^ME[y, the 
slopes for various S are similar, the efficiency can be described by 
e PC This is flatter than previous estimations of e oc M in 
iNaravan. Mahadevan & Ouataert jl998h where 6 = 10"'. This dis- 
crepancy is not because of 6 or outflow effect, but seems to be sim- 
ply becaus e that the estimation in previous work is "rough". From 
Figure 7 in lNaravan. Mahadevan &Ouataertl ( ll998l) . as the accre- 
tion rate changes by 3 orders of magnitude, i.e., from 10"'*MEdd to 
10"'MEdd> the bolometric luminosity varies by 5 orders of mag- 
nitude, i.e., from 10"'LEdd to 10"^LEdd- So we should also have 
6 oc m" *""-^ , fully consistent with our result. 

• In the ADAF regime, the radiative efficiency strongly depends 
on the value of S. This is because a larger d implies that more energy 
will be received by the electrons, subsequently higher radiative ef- 
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ficiency. But note that when M is small, the efficiency is still very 
low. In the case of the accretion flow in Sgr A*, if we adopt the def- 
inition of Eq. ([8}, the radiative efficiency will be ~ 4 x 10"^, which 
is lower than that of a standard t hin disk by a factor of 4 x 10""* 
jYuan. Ouataert & Naravaij|2003h. Since ^p^, = i^son d ~ 1 0^ 
and ~ 0.3 are adopted in I Yuan. Ouataert & Naravai] 1 2003h . the 
mass loss in outflow contributes (1/10^)""^ ~ 0.04, the other factor 
(4 X 10"''/0.04 a 10"-) is because of energy advection by both ions 
and electrons (refer to section |4j2}. 

• When M„et ~ Mcadaf^ for different 6 the slopes are all very 
steep and the values of e become comparable. This is because in 
this regime of M, ^vis.e is compensated by ^ie(~ i?vis,e) in the elec- 
tron energy equation. Moreover, the main radiative process is the 
Comptonization of synchrotron photons. This process is very sen- 
sitive to the optical depth, or accretion rate, of the accretion flow. 
This is why we have a steep slope. 

• When 6 = 0.5, the radiative efficiency can be as high as 3%, 
even when M|,et is as low as ~ 2 x lO^^Medd- Of course, if the 
definition of Eq. {Sj is adopted, the efficiency will be lower by a 
factor of (Rs/RoutY (~ 0.16 in our case). Still, this implies that the 
efficiency of ADAFs is not as low as people sometimes imagine. 
For other values of i5, when M„et ^ Mc,_adaf, i.e., in LHAFs, the ra- 
diative efficiency is all quite high. One application of this result is 
that we should not observe any large luminosity change during the 
state transitions from hard to soft in black hole X-ra y binaries. This 
is we ll consistent with X-ray observations (e.g., IZdziarski et al.l 
l2004h . 

• The radiative efficiency of the two-phase accretion model is 
nearly independent of the accretion rate, similar to the standard thin 
disk model. The radiative efficiency can be as high as ~ 8% in our 
definition. This is close to or even slightly larger than that of the 
standard thin disk (note we consider a Schwarzschild black hole). 
For two-phase accretion flows, the electron advection is zero. Sum- 
ming up equations ^ and (O, we have q^^i = ^vis - ^adv.i. Since 
both ^vis and q.^^^j are proportional to density thus M, e will be a 
constant (also independent of S). As q^^ii^j < for a LHAF while 
lidv.i = for a SSD, the efficiency of LHAF can be slightly higher 
than that of a standard thin disk, if the outflow effect is not consid- 
ered. 

• The range of accretion rate within which two-phase accretion 
model exists spans a factor of 8 - 10 for our chosen parameters. 
This range is roughly independent of 6. 

For the convenience of use, in the following we provide a 
piecewise power-law fitting to efficiency for ADAF and Type I 
LHAF. We assume. 



e(M„et) = eo 



(11) 



where the normalization = lO^^Medd = 10"'LEdd/c^. The fit- 
ting results can be found in Table [T] The boundary accretion rates 
are adjusted after the fitting so that a continuous fitting function is 
achieved. 

So far our numerical calculations and fittings are only for a = 
0.1. In the literature, o- = 0.3 is also widely adopted. We therefore 
have also calculated the a = 0.3 and S = 0.5,0.1 and 10"' cases. 
We find that the following formula, with coefficients taken from 
corresponding o- = 0. 1 cases, presents good fit to the cases of d = 
0.5 and 0.1, 



Table 1. Piecewise power-law fit formulae of radiative effi- 
ciencies for ADAF and Type I LHAF. 



Cases 


M„et/MEdd Range 




eo 


index a 


6 = 


0.5 


<2.9x 10-' 




f.58 


0.65 






2.9 X 10"^ - 3.3 X 10" 


-3 


0.055 


0.076 






3.3 X 10"' -5.3x 10" 


-3 


0.17 


1.12 


S = 


0.1 


<9.4x fO-' 




0.12 


0.59 






9.4 X fO-5 -5.0x fO- 


-3 


0.026 


0.27 






5.0 X fO-' - 6.6 X fO" 


-3 


0.50 


4.53 


S = 


10-^ 


<1.6x IQ-' 




0.069 


0.69 






1.6 X 10"' -5.3x 10" 


-3 


0.027 


0.54 






5.3 x 10"' -7.1 X 10" 


-3 


0.42 


4.85 


5 = 




<7.6x 10"' 




0.065 


0.71 






7.6 X 10"' - 4.5 X 10" 


-3 


0.020 


0.47 






4.5 X fO"' -7.f X fO" 


-3 


0.26 


3.67 



Note: See context for the definition of radiative efficiency e. The fitting 
takes the form e(M„ct) = eo (^net/^r ) ■ where the normalization accretion 
rate is fixed at = IO^^Mejj = lO^'LEjd/c-. For other values of a, the 
fitting formulae can be found in Eq. )121 or Eq. U3t . depending on S and 
accretion rate. 



+ 

S 

"> 
a- 



III 



1.0x10:^ 




1.0x10-3 



10 100 

R/R. 

Figure 2. The advection factor of electrons (/e) when Mq < 10"'MEdd for 
6 = 0.5. The accretion rate at fiom = iO~R^ for each curve is labeled in the 
plot (in unit of MEdd)- 

here = 0. lo-MEdd = Q'^Edd/f"- Note that for o- = 0.1, it recovers 
the previous definition. For d = 10"' case, if M„et ^ 7 x 10"-Q'^MEdd 
(the typical ADAF regime), the efficiency can be nicely fitted by 
the following formula, 



e(M„e,) = eo 



(13) 



When the accretion rate is higher, especially when M„et ^ 5 x 
10"2Q'MEdd (transition from ADAF to Type I LHAF), Eq. ^ 
presents good fit to the data. 

4.2 Energy balance for ions and electrons 

We now investigate the energy balance relationship among the 
terms in the energy equations of ions and electrons. We define the 
advection factors for ions and electrons as follows. 



(12) 



,/e 



?adv,i _ 



1 



(1 - S)q,, 



1 



9rad 



(14) 

(15) 
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Note that from Eq. l ll4b . one can get, 



?1C 



(1-^X1-/0 



(16) 



which characterizes the relative importance of i^vis.c and in the 
energy equation of electrons. Evidently from this equation, we find 
that for large value of 5 (e.g., S ^ 0.1), viscous heating to electrons 
will be the main heating term to the electrons, provided that the ions 
are advection-dominated (i.e., /i > 0.9). Even for lower 6, ^vis.e will 
still be the main heating term at low accretion rate, where /i = 1. 

Below we focus only on the case of ^ = 0.5 (the long-dashed 
curve in Fig.[T]l, since this value of S is most favored theoretically 
from the detailed study of Sgr A* jYuan. Ouataert & NaravanI 
l2003h . The results are shown in Figs.[2]&[3] Note that q^i^.i = Ins.s 
for this choice of 5. We describe the results at the following three 
regimes of accretion rate. 



• Mo<3.0x 10-^ Mem. 
In this regime, both the ions and electrons are strongly advection 
dominated, i.e., = 1,/e = 1. Both the Coulomb coupling and 
the radiative cooling rate are negligible compared to the viscous 
heating rate. Specifically, we have. 



ions : ^adv.i ~ gv,s,i » i?ic, 
electrons : ^^jv.c ~ ^vis.c » ?rad&9ic- 



III 




Figure 3. The energy balance relationship for ions and electrons when 
Mo>5 X lO^^M^ad for 6 = 0.5. Upper panel: advection factor of ions 
//. Lower panel: advection factor of electrons /e (black) and ^vis.c/^ie (red). 
The lines with the same type have the same accretion rate with the upper 
panel. The curve of ^vis.c/'Jie for ^0 = 5.0 X W^^Meu is not shown here 
since its value is too large. 



(17) 



• 3.0 X 10-^MEdd < Mo < 1-0 X lO-^Medd- 
In this regime, the accretion rate is still low enough, the ions re- 
mains advection-dominated (/i « 1). For the electrons, the radia- 
tive cooling rate is high because of the high density of the accre- 
tion flow. As illustrated in Figs.[2]&[3] the electrons are radiation- 
dominated, i.e. the advection factor /e ~ in the inner R < 20Rs 
regions. 

ions : g^dv., ~ gv.s.i » gic, 
electrons : q„i a; ^^is.c » ?adv.c&9ic- (18) 

. 1.0 X 10-2 MEdd <Mo< 3.2 X 10-2 MEdd- 
This is the transition regime from an ADAF to a Type I LHAF. As 
the accretion rate (or the density) increases, the Coulomb coupling 
between ions and electrons becomes so strong that the ions are no 
longer advection-dominated, i.e. q,c ~ q^isj- Above Mc^adaf, we 
have > ^vis.i in some regions of the hot flow, which means that 
the flow enters into the Type I LHAF regime. All the energy terms 
in the ion energy equation are roughly comparable to each other. 
The electrons radiate away nearly all the energy they receive via 
viscous heating and Coulomb collision, i.e. = 0. The relative 
importance of ^vis.e and q[c depends on various parameters, i.e. Mo 
and i5, and the radius R. 

ions : g^jv,. ~ gv.s.i ~ g,e, 
electrons : q,^^i a; (^v.s.e + ?.c) » ?adv,c- (19) 

For other choices of a, the above arguments still hold, except 
that the accretion rate regime is replaced by (a/O.l) times the value 
listed above. 

Above Mcr,LHAF. purely hot solutions do not exist; and the ac- 
cretion flow enters the two-phase regime. In this regime, we expect 
that /e - 0, i.e. q,^ = ^vis.e + qm- 



5 DISCUSSIONS: CAVEATS AND RADIO-X-RAY 
CORRELATION 

In our calculations, we only consider the local Compton scat- 
tering, namely the scattering between photons and electrons oc- 
curred at the same region where the photons are produced. 
However, since a hot accretion flow is usually optically thin 
in the radial direction, the photons produced at one certain ra- 
dius can in principle travel for a long distance and collide with 
electrons at another radius. Such a "global" Compton scatter- 
ing effect has been systematically investigated in previous works 
llPark& Ostriker2001, 2007:'Yuan. Xie & Ostrikerl2009l : IXie et al.l 
[2010; Niedzwiecki, Xie & Zdziarski 2012) . It was found that it 
plays a significant cooling and heating roles in the region of 
50R,<R<W0R, and ie>5 X 10^^?^, respectively, when the ac- 
cretion rate is high enough so that t he total luminosity emitted 
from the accret ion flow Lboi ^ 2%Lem jYuan. Xie & Ostrikeil2009l : 
IXie et alj|2O10l) . One consequence is that the radiative efficiency 
will be lower by a factor of 2 compared t o the case that this effect 
is not taken into account jXie et al1l2010l) . In addition, the highest 
luminosity a hot accretion flow can emit will be constrained to be 
Lboi ^ l^^Edd- Above this limit, the global Compton cooling and 
heating will be so strong that n o steady hot solution can be found 
and the system will "oscillate" jYuan. Xie & Ostrike3l2009l) . Note 
that if Ltioi < 2%LEddi or outer boundary radius of hot flow is small, 
i.e. i^oui ^ 50 - 100i?s, the global Compton scattering effects will be 
unimportant. The latter is the case of luminous hard state of black 
hole X-ray binaries. 

In our model, we assume that the magnetic filed is tangled and 
weak, thus it does not play any dynamical role. Numerical simula- 
tions have shown that a large-scale toroidal magnetic field is likely 
to exist in the inner regio n of the accretion flow, imposed on the 
stochastic component (e.g. lHirose et al.ir2004l) . T he effect of such a 
field has been studied by self-similar approaches jAkizuki & Fukuel 
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bOOdlAbbassi, Ghanbari & Naiiaill2008l;lBu, Yuan & Xiell2009l) or REFERENCES 



global calculations i Oda et al] 20071 20121) . Especially, the global 
solution with strong large-scale magnetic fields indicates an in- 
crease in the hig hest luminosity a hot accretion flow can achieve 
dOda et alj20Ta) . 

Throughout this paper, we fix the outer boundary condition 
(T\, Tc, v) in our calculations. The effect of outer boundary con- 
dition on t he dy namics of accretion flow has been studied in 
lYuan et al] ( I2OO0I) . Obviously, it will also influence the radiative 
efficiency. 

A correlation between the radio and the 2-10 keV luminos- 
ity (Lx) has been found among the hard state of black hol e X- 
ray binaries a nd low-luminosity active galact i c nuclei (Corbel et al.l 
200d. l2003l: lOallo. Fender & PoolevI 120031; iMerloni et al.i |2003| ; 



Gallo. Miller & Fendeil2012l : ICorbel et al.l2012h . which is well de- 
scribe d by a power-law, LRj dio «: with index p ~ 0.6. Re- 
cently IZdziarski et alj ( 12011 ) found that this correlation extends 
to intermediate and soft states for Cyg X-1, if only the luminos- 
ity from hot disc is used. This correlation ha s been quantitativel y 
explained by the coupled jet-ADAF model in I Yuan & Cuil ( l2005h . 
in which the radio and X-ray emissions are dominated by the ra- 
diation from the jet and ADAF, respectively. It is interesting to 
note that there are now growing number of sources which show 
that when Lx > 4 x 10^* erg s"' the radio/X-ray correl ation follows 
a steeper power-law, with index p ~ 0.98 or 1.4 jCoriat et al.l 



l201lHGallo. Miller & Fendeil2012feo"rbel et al.l2012t) . Below this 
critical luminosity, the sources return to the ~ 0.6 correlation at 
~ 10" ergs"'. Between 4 x 10'* ergs"' and 10^^ erg s~', th e ra- 
dio luminosity remains almost unchanged. ICoriat et"al] ( 1201 ih pro- 
posed that one way to explain the steep (p ~ 1.4) correlation is that 
the radiative efficiency of the hot accretion flow is independent of 
the accretion rate, if the ratio of the mass loss rate in the jet and the 
accretion rate in the accretion flow is a constant. As shown by Fig. 
1, this is the case for our two-phase accretion flow (Type II LHAF). 
Moreover, as also shown by this figure, the efficiency curve of Type 
I LHAF is very steep, which means that a small change of accretion 
rate will result in a large change of Lx. This feature is obviously at- 
tractive to explain the "flat transition" between 4 x 10'* erg s"' and 
10'^ ergs"'. The reason why some sources follow a single ~ 0.6 
correlation while others follow three branches is simply because of 
different values of a among these sources. If a is large, adaf 
will be large, so no transition to LHAFs will occur throughout the 
evolution of M during the outburst. This is why we can only ob- 
serve one single ~ 0.6 correlation. If on the other hand o- in a 
source is small, Mcf.adaf will be small thus the sources will en- 
ter the two-phase LHAFs regime during the outburst. In this case, 
three branches of correlation should be expected. In a future work 
we plan to investigate the correlation in detail. 
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